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Abstract 



We construct consistent theories of multiple interacting spin-2 fields in arbitrary 
spacetime dimensions using a vielbein formulation. We show that these theories have 

_ 1 the additional primary constraints needed to eliminate potential ghosts, to all orders in 

the fields, and to all orders beyond any decoupling limit. We postulate that the number 
of spin-2 fields interacting at a single vertex is limited by the number of spacetime 

l/~) dimensions. We then show that, for the case of two spin-2 fields, the vielbein theory is 

equivalent to the recently proposed theories of ghost-free massive gravity and bi-metric 

O 

£Sj gravity. The vielbein formulation greatly simplifies the proof that these theories have 



X 
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an extra primary constraint which eliminates the Boulware-Deser ghost. 
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1 Introduction and summary 

Theories that propagate a massive spin-2 particle have traditionally been plagued by the 
presence of an additional unstable degree of freedom - the Boulware-Deser ghost [1]. For 
the case of a single interacting massive graviton, this problem was solved only recently, by 
the dRGT (de Rham, Gabadadze, Tolley) massive gravity theories described in [2, 3]. These 
theories were arrived at by choosing interaction terms so as to raise the cutoff of the effective 
theory [4, 5], after which the Boulware-Deser ghost is automatically vanquished. (See [6] for 
a recent review of these aspects of massive gravity.) 

The dRGT theories can be put in the following form [7] , 

2 4 " 



S = — — / d^x \f—g 
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n=0 



(1.1) 



This is the Einstein-Hilbert kinetic term for the metric g^, supplemented by a potential 
term containing no derivatives of the dynamical metric g^ but depending explicitly on an 
external flat metric rj^ which breaks the diffeomorphism invariance of the m = theory 1 . 
Here, S n is the n-th elementary symmetric polynomial (see Appendix A) of the eigenvalues 
of the matrix square root of g^rj^, and the j3 n are free coefficients. 

The dRGT theories (1.1) are ghost-free [8, 9, 10, 11, 12, 13, 14, 15]. This means they 
propagate only the five degrees of freedom of a massive spin-2 field, and not the additional 
unstable degrees of freedom which plague all interaction terms not of the form of those in 
(1.1) (including those studied in [1]). 

The counting of degrees of freedom can be done in the Hamiltonian formulation. As 
in General Relativity (GR), we adopt ADM variables [16] and Legendre transform with 
respect to the spatial components of the metric, obtaining a 12 dimensional phase space 
consisting of the 6 components of the symmetric spatial metric g^, and the 6 components 
of its canonical momenta n^ . In massless GR, the lapse and shift variables appear as 
Lagrange multipliers, enforcing 4 first class constants which remove 8 degrees of freedom, 
leaving a 4 dimensional physical phase space which describes the 2 polarizations of the 
massless graviton and their conjugate momenta. For generic potentials however, the lapse 



1 Of course, the diffeomorphism invariance of massive gravity can be restored by the introduction of 
Stiickelberg fields n . In this work we treat massive gravity in the unitary gauge in which 7r = 0. 



and shift variables appear as auxiliary variables rather than Lagrange multipliers. They 
can therefore be eliminated by their own equations of motion, leaving no constraints, so the 
entire 12 dimensional phase space is physical, describing the 5 polarizations of the massive 
graviton and the single Boulware-Deser ghost, as well as their conjugate momenta. For the 
special interactions of dRGT massive gravity, however, the lapse still survives as a Lagrange 
multiplier, enforcing a single primary constraint [8]. This leads to a secondary constraint 
[11] which, together with the primary constraint, forms a second class pair of constraints. 
This leaves a 10 dimensional physical space space, just right for the 5 degrees of freedom of 
a massive graviton and its canonical momenta. 

The potential of (1.1) is awkward, involving matrix square roots. In particular, this 
makes it difficult to identify the extra primary constraint which removes the Boulware-Deser 
ghost. In the formulation (1.1), this constraint is only seen after a complicated re-definition 
of the shift variable. The square root structure suggests that vierbein variables E A = E^dx 11 
may be better suited for describing this theory since they are, in a sense, like the square root 
of the metric, 

g lw = E li A E v B riAB. (1-2) 

In this paper, we will show that the dRGT theory (1.1) is precisely equivalent to the following 
action written in terms of vierbeins, 

m 2 ( r 

S = -^-l d 4 x (det E)R[E] 
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where 1 A = 6 A dx^ is the identity vierbein, which can be thought of as a vierbein for the 
flat background metric r/^. The flat space epsilon symbol is denoted by €a 1 a 2 —Aa' 

The potential terms now appear simply as all possible wedge products of the vierbein 
of the dynamical metric with the vierbein for the background. The structure of the wedge 
product is responsible for the appearance of the symmetric polynomials in (1.1). In these 
variables, we will find that it is almost trivial to see the existence of the extra primary 
constraint which makes the theory ghost-free, to all orders in the fields, and to all orders 
beyond any decoupling limit. 

In (1.1), the flat metric rj^ is fixed. One can promote this to a general reference metric, 



ffj, v , and the theory is still ghost-free [17]. Going further, one can promote / M „ to a dynamical 
metric by adding a kinetic term for it, 



S = (fx 



Ml ^_ Ml ,— m 2 M 2 ^* , 
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n=0 

Unlike the theory with a fixed reference metric, there is now an overall diffeomorphism 
invariance. This is a bi-gravity theory, describing at linear level a massless helicity-2 graviton 
plus a massive spin-2 graviton, for a total of 7 degrees of freedom. For a general choice of 
interaction terms, a generic bi-gravity theory would also have a Boulware-Deser-like ghost, 
for a total of 8 degrees of freedom. But with the interaction terms of (1.1), it can be shown 
that the ghost is absent [18] so that there are in fact 7 degrees of freedom. 

We will show in this paper that the bi-metric theory can also be written simply in 
vierbein form, 

Ml f , M 2 f . 

S = — 9 - / d*x (det E w ) R[E W ] + — j - / d 4 x (det E m ) R[Epl 

m 2 M 2 f fl 
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where E&) and E( 2 ) are vierbeins for the two metrics g^ v and f^ u , respectively. In this form, 
the extra primary constraint responsible for eliminating the Boulware-Deser ghost will be 
easy to see. 

By using the vierbein formulation, we can go beyond the bi-gravity theory and construct 
ghost-free interaction terms that directly mix more than two gravitons. These are terms 
which cannot be easily inferred using the metric formulation. We can have tri-metric vertices 
which link together 3 different metrics. There are three possibilities, depending on which 
vierbein appears twice in the wedge product, 

£a 1 a 2 a 3 a 4 E(i) 1 A E(i) 2 A E( 2 ) 3 A E(3> 4 , 
£a 1 a 2 a 3 A4, E(i) x A E( 2 ) 2 A E( 2 ) 3 A E( 3 ) 4 , 
~e Al A 2 A 3 A 4 E (1) Al A E (2) A2 A E (3) A3 A E (3) A4 . (1.6) 

Finally, there is one possible tetra-metric vertex which can mix 4 different vierbeins, 

e AlA2AsA4 E w Al A E (2) A2 A E (3 ) A3 A E w Ai . (1.7) 



There is no ghost-free vertex of this type which can link 5 or more metrics, since this would 
require wedging 5 or more vierbeins, which can't be done in 4 dimensions. 

Generic theories of M interacting spin-2 fields will contain the degrees of freedom of 1 
massless spin-2, J\f — 1 massive spin-2's and M — 1 Boulware-Deser-like scalar ghosts. This 
is because each metric carries 6 potentially propagating degrees of freedom, rather than 5, 
and there is only one overall diffeomorphism invariance. In this paper we will show that the 
above potential terms have the additional primary constraints needed to eliminate the M — 1 
ghosts, non- linearly to all orders in the fields and to all orders beyond any decoupling limit. 

The vierbein formulation and the proof of ghost-freedom can be extended to D space- 
time dimensions. We will thus work in arbitrary dimensions for the remainder of the paper. 
The dRGT interaction terms remain the same, but with D symmetric polynomials. For a 
theory with multiple interacting gravitons indexed by /, our conjecture is that the following 
term is the most general ghost-free potential in D spacetime dimensions, 

AT 
U= Yl T hh - lD e AlA2 ... AD E (h) M A E ( z 2 ) M A . . . A E (Id) Ad , (1.8) 

h,...,i D =i 

where T IiI2 '" Id is a completely symmetric constant tensor of coefficients. We will show that 
this term yields the primary constraints needed to eliminate all of the M — 1 ghosts, in all 
dimensions, and to all orders in the fields and beyond any decoupling limit. Since these 
terms are constructed by wedging together combinations of the vielbeins, they are limited 
by the number of spacetime dimensions. These ghost-free vertices can directly connect only 
up to D different gravitons. 

Following the ideas of [4, 19], we can represent these multi-metric theories graphically. 
The bi-metric theory is given in figure 1. The two nodes represent the two metrics, each of 
which comes with an Einstein-Hilbert term. The line connecting the nodes represents the 
bi-metric mass term, which mixes the two metrics. In the vielbein formalism, each node has 
a separate diffeomorphism invariance, and a separate local Lorentz invariance (LLI). The 
mass term breaks this down to the diagonal subgroup of diffeomorphisms and LLPs, those 
transformations which act the same way on both metrics together. 

Continuing in this way, we can link up any number of gravitons using bi-metric inter- 
actions, for example as in figure 2. There will be one unbroken diffeomorphism and one LLI 
for each "island" (i.e., disconnected subgraph) in the graph. The spectrum will consist of 
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Figure 1: Visual depiction of a bi-metric theory. The two circles represent the two vielbeins, the 
bar between them the interaction terms which mix them. 

one massless graviton for each island, and the rest of the gravitons will be massive. (This 
does not contradict theorems forbidding multiple interacting massless gravitons [20], since 
the islands don't interact with each other.) We will show that an arbitrary graph, when the 
links are constructed using the bi-metric interactions of (1.5), has the primary constraints 
to make it ghost-free to all orders. Thus it is now possible to build the kinds of theories in 
[4] , for deconstructing gravitational dimensions, in a ghost-free manner. 




Figure 2: Visual depiction of a general multi-metric theory with only bi-metric interactions. 

The new multi-metric interactions discussed in this paper can be represented as nodes 
where three or four metrics meet, as in figure 3. These can now be used as ingredients in a 
general theory graph, and the result will still be ghost-free. 





Figure 3: Visual depiction of the new multi-metric ghost-free interactions discussed in this paper. 



Vielbeins have proven useful in the study of massive gravity. The vielbein formalism 
within the context of massive gravity has been developed in [21, 22, 23, 24, 12]. Here we 
suggest that they may very well be essential. As we'll see, there is no straightforward way 
to write down an equivalent metric formulation for theories including these tri-metric or 
tetra-metric vertices, or even theories with only bi-metric interactions whose theory graphs 
contain closed loops. The vielbein theories provide a new and very wide parameter space 
in which to search for theories that may address issues such as superluminality [25] and 
the low cutoff of the effective field theory [4]. There are parallels in galileon theory which 
indicate this might be possible. For example, single field galileon theories [26] always have 
superluminality [26, 27], whereas in multi-field galileon theories [28, 29], the parameters can 
be chosen so that there is no superluminality [30]. 

In this paper our analysis will proceed from the general to the specific. We first intro- 
duce the general interaction terms in vielbein form and perform the Hamiltonian analysis. 
We show that the vielbein formulation makes it straightforward to identity the primary 
constraints which eliminate the potential ghosts. We then consider the recently proposed 
ghost-free bi-gravity and massive gravity theories which were formulated in terms of metrics 
rather than vielbeins. We show that these theories are dynamically equivalent to the vielbein 
multi-gravity theories, for the case of two metrics. We repeat the Hamiltonian analysis for 
the case of massive gravity in the vielbein formulation, addressing the subtleties for this case. 

Conventions: D refers to the number of spacetime dimensions while d = D — 1 is the 
number of spatial dimensions. We use Greek letters /x, u, etc., for spacetime indices and i, j, 
etc., to denote their spatial components. These are raised and lowered with the full metric, 
g^ v and the spatial metric g^ respectively. We use capital letters A, B, etc., for Lorentz 
indices and lower case letters a, b, etc., for their spatial components. These are raised and 
lowered with the full Minkowski metric, i]ab and the flat spatial metric 8 a b, respectively. The 
indices (J), (J), etc., label the M spin-2 fields. The Einstein summation convention is used 
for all indices but the (I), (J), etc. We use the mostly plus metric signature convention, 
Vij.v — ( — )+)+;+)■■■)■ E denotes the D-dimensional vielbein while e denotes the spatial 
vielbein. Indices are (anti-) symmetrized with weight 1, i.e. [/iz/] = |(/iz/ — vfi). The epsilon 
symbol has a tilde, and is always defined so that I012-D = +1, regardless of the placement 
of indices or the signature of the metric (so i 012 '" D = +1 as well). 



2 Mult i- metric theory 

We start with the general case of multiple interacting spin-2 fields in D dimensions. We have 
a collection of M metrics gw^, indexed by /, each with its own Einstein-Hilbert term and 
Planck mass Miry, 

yr D-2 

C = J2^^V^)Ri9(n}. (2.1) 
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Without interaction terms, this action has M separate diffeomorphism invariances, one for 
each /, given by arbitrary functions fiif(x), 

avw -> ~dxy~dx^ 9(I)af} (/(/)(x)) • ( > 

As we wish to express the interactions in vielbein form, for each metric 9(i)u, u (x) we 
introduce a vielbein one-form E</) = E(i) (x)dx fl , which relates the metric to the flat 
Minkowski metric 77 ab, 

«V = %V^' ( 2 - 3 ) 

We denote the inverse vielbein by E(rf t A (x), so that Eii^Eii)^ = 8 A and Eiif" A E(i) v = 8%. 

Without interaction terms, we now have, in addition to the diffeomorphism invariances, 
a local Lorentz symmetry for each J, 

E(% ->■ A(/) B E ( i)^ , (2.4) 

where A(j) B {x) G SO(l,D — 1). The inverse transforms as E(if A — > A(i) A Eiif B , where we 
define the inverse At/j^j 8 = (A(i) -1 ) so that A(z) c A(j) B = A(i) cAfi)^ 7 = 5 B . The vielbeins 
transform as one-forms under the local diffeomorphism symmetries, 

Ei^W^^E^iMx)). (2.5) 

Note that, in pure GR, the passage from metric to vielbein via the replacement (2.3) is 
nothing but the Stiikelberg trick: the metric is symmetric and so has D(D + l)/2 components, 
while the vielbein has no restriction other than invertibility and so has D 2 components. We 
have thus introduced D(D — l)/2 new, unphysical gauge invariances (the local SO(l, D — 1) 
symmetries) along with D(D — l)/2 new unphysical fields. 



2.1 The interaction terms 

Our goal is to write down the most general interactions among these spin-2 fields which do not 
introduce ghosts. We look for interactions which would be most relevant at long distances, 
those involving no derivatives of the metric or vielbein. In addition, we demand that at least 
one overall diffeomorphism and one overall local Lorentz symmetry remain unbroken, that 
is, the one for which the f(ij l {x) are equal for all J, and the A(/) A ,(x) are equal for all /. 
This requires the potential to be a diffeomorphism and Lorentz scalar constructed out of the 
vielbeins. 

Our conjecture is that the following term is the most general ghost-free potential of M 

spin-2 fields, 

N 
u= J- T hh - lD i Al A 2 -A D Edo Al A E(, 2) M A ... A E (lD) A ° , (2.6) 

h,...,i D =i 

where T IiI2 " Id is a completely symmetric constant tensor of coefficients, and £a 1 a 2 -~a d is 

the totally anti-symmetric flat space epsilon symbol. This term is invariant under an overall 

local Lorentz transformation and transforms as a scalar density of weight 1 under general 

coordinate transformations. In D dimensions, the total number of different terms in this 

( M + D-l \ 

A/ -metric interaction is ] , the number of independent components of a rank 

D symmetric tensor in M dimensions. 

Let's look at some of the possible terms. The simplest case is when there is only one 
metric g^ u , and hence one vielbein E A . Then there is only one possible term, 

e AlA2 ... Ao E Al A E A2 A . . . A E A ° , (2.7) 

corresponding to a cosmological constant, oc detE oc y/—g. 

If there are two metrics gw and (?(2>, in addition to a cosmological term for each metric 
separately, there are now D — 1 possible terms which mix the two, 

e Al A 2 A 3 -A D E( 2 ) Al A E (2) A2 A E (2) Aa A ... A E (2) Ad , 
e Al A 2 A 3 -A D E (1 ) Al A E (2) A2 A Et/ 3 A ... A E (2) A ° , 
e Al A 2 A 3 -A D E(d Ai A E W A2 A E (2) Aa A ... A E (2) A ° , 

e Al A 2 A 3 -A D E(d Ai A E (1) fl2 A E w As A ... A E W A ° . (2.8) 



We will see that these are equivalent to the interaction terms of ghost-free bi-gravity, and 
when one of the metrics is non-dynamical, ghost-free massive gravity 

With multiple metrics, we are allowed to wedge together all the various metrics in all 
possible ways. Let us enumerate all the possibilities for D = 4. With one metric, there is 
only the cosmological constant term. For two metrics, there are the 5 terms listed in (2.8); 
the two cosmological terms for each metric separately plus the three mixing terms. For three 
metrics, there are 3 cosmological terms, 3 terms mixing each of the 3 possible pairs of metrics 
(for a total of 9 bi-metric terms) and there are 3 terms that mix together all three of the 
metrics, 

£a 1 a 2 a 3 a 4 E(i) * A E(i) 2 A E( 2 ) 3 A E(3> 4 , 
£a 1 a 2 a 3 a 4 E(i) x A E( 2 ) 2 A E( 2 ) 3 A E( 3 ) 4 , 
~e AlA2 A 3 A 4 E (1) Al A E (2) A2 A E (3) A3 A E (3) A4 , (2.9) 

for a grand total of 15 terms. For 4 metrics, there are 4 cosmological terms, 3 terms mixing 
each of the 6 possible pairs of metrics (for a total of 18 bi-metric terms), 3 terms that mix 
together all three metrics in each of the four 3-metric subsets (for a total of 12 tri-metric 
terms), and finally one term which mixes all four metrics, 

e AlAaAaAi E m Al A E (2) A2 A E( 3) A3 A E w Ai , (2.10) 

for a grand total of 35 terms. With 5 or more metrics, there are no terms which mix 
together all of the metrics, so there are only the terms which mix 4 or fewer, for each 

subset of 4. In general, the total number of terms in an A^-metric theory in 4 dimensions, 

/ 7V + 4-1 \ 
is ] , which is the number of independent components of a rank 4 symmetric 

tensor in TV dimensions. 

For general dimension D, with M metrics, the number of possible terms which mix 
together n of the metrics (n < A/") can be written 

(2.11) 

The first factor is the number of ways of choosing the subset of n metrics which are coupled 
together, and the second factor is the number of ways of partitioning the D terms in the wedge 

9 




/ Af + D- 1 
product among the n metrics. Summing over n — 1 • • ■ N , we find a total of 

possible terms. Note that if J\f > D, there are no terms which mix together all the metrics. 

We can graphically represent a theory with TV-metrics and the interaction term (2.6) 
by drawing a "theory graph," (or gravitational quiver diagram) following [4, 19, 31, 32]. For 
each metric, we draw a dot. For each term in the interaction potential (2.6) involving only 
two of the metrics, we draw a line connecting the corresponding two dots. For the terms 
involving three or more metrics, we draw a cubic, quartic, etc. vertex as in Figure 3, and 
connect its edges to the corresponding metrics. The number of unbroken gauge and local 
Lorentz invariances is the number of disconnected "islands" in the theory graph. 

It will sometimes be more convenient to work with a matrix expression for the the 
interaction terms, rather than the wedge products. Accordingly, we derive the following 
useful result 

£a 1 a 2 -a d ~E'(i 1 ) 1 A E(i 2 ) 2 A ... A E(/ d ) d 



D 



= (det E (h) d D x) i Al A 2 -A D e B ^- B ° 6£ (E^E^)^ ■ ■ ■ (E (Ii r l E (lD) ) B A 

= (det E (h) d D x) 5(1, E (Il f l E V2h ..., E (h) - l E {lD) ), (2.12) 

where (E^) -1 E(i 2 )) B 2 = E^f^ E(i 2 ) 2 and dPx = dx 1 Adx 2 A. . .Adx D is the volume element, 
and S is the multi-matrix symmetric polynomial defined and discussed in Appendix A. 2. 

To see this, start by writing the wedge products in terms of epsilon symbols and the 
volume element, 

~e AlA2 ... AD E (h) M A Eu/ 2 A ... A E (Id) Ad = 1 Ai a 2 -a d l^~^ E ih)fI Al E^ ■ ■ ■ E (lD)/ ff d D x . 

(2.13) 
Then single out E(i x ) and re-express the curved space epsilon symbol in terms of the flat 
space epsilon symbol, the determinant detEii^ and a bunch of inverse vielbeins EnJ 1 ^ 

(det E (h) d D x) € AiA2 ... Ad ~e B ^- B ° E (h) ^ E ih) £ ■ ■ ■ E {Id) ^ E {h) % E w % ■ ■ ■ Bk h f§ D 
= (det E ih) d D x) i Al A 2 -A D ~e B ^~ BD 8 B Al {E (H) - l E (l2) ) A2 ■ ■ ■ (E^E^)^ . (2.14) 

This expression singles out one of the vielbeins to be in the determinant, however, this choice 
is arbitrary - there are other equivalent expressions depending on which vielbein is chosen 
to be in the determinant. 

10 



As an example, consider the interaction terms in D = 4. Define the matrices 

X = E(i) E(2), Y = E(i) E(3), Z = E(i) E(a). (2.15) 

The ghost-free potentials take the form 

l MMAzM E W M A E W M A E W A3 A E (2) A4 = 6 (det E (h) d D x) [X] , 
i Al A 2 A 3 A 4 E (1) Al A E W A2 A E (2) A:i A E (2) A4 = 2 (det E (h) d D x) ([X] 2 - [X 2 ]) , 

£ AlA2A3 A 4 E (1) Al AE W A2 AE (2) As AE (3) A * = 2 (det E (h) d D x) {[X][Y] - [XY]) , 

l AlMA3M E w A ^ AE {2) M AE^ AE {2) M = (det E ih) d D x) ([X] 3 - 3[X][X 2 ] + 2[X 2 

e^^^^Et^AE^AE^AE,^ = (det E (lx) d D x) ([X][Y][Z] - [X][Y Z] 

-[Y][XZ] - [Z][XY] + [XYZ] + [XZY] 

(2.16) 

as well as every non-redundant permutation of e(i), ep), e<3) and e<4> Along with the four 
cosmological constants, these give the 35 interaction terms described in the previous section. 

2.2 Hamiltonian formulation 

In this section we perform a Hamiltonian analysis of the multi-vielbein theory with a general 
interaction term (2.6). To this end, we perform a d+ 1 decomposition of the vielbein into 
canonically conjugate ADM variables. 

A general vielbein can always, by a local Lorentz transformation, be put into upper 
triangular form (upper triangular vielbeins will be written with a hat), 



N N%« \ ^ I jj 

'** In „ a ' A \ N 



Here the X and X* are the .D time-like components. The spatial vielbeins e" 1 contain (D — l) 2 
components and are related to the spatial part of the metric by g^ = e i °e J - 6 5 6. By writing 
out the metric of this vielbein, we see that X and X 1 are the usual lapse and shift of the 
metric ADM decomposition [16], 

a - R I -N 2 + N'Ni X . 

g, u = E A E u B r ]AB =-[ ) . (2.18) 

Nj Qij 
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The upper triangular form does not completely fix the local Lorentz invariance. It 
leaves a residual local spatial rotation. There are D components in the N, N l and (D — I) 2 
in the spatial vielbein. The remaining D — 1 components of the general vielbein have been 
fixed by using the upper triangular gauge choice. 

We can formulate an arbitrary vielbein as the action of some standard boost on an 
upper triangular vielbein 2 . For every given d- vector p a , we define a standard Lorentz boost 

Pb K + -tfjfpb J 

where indices on p a are raised and lowered with 5 a b and 

7 = y/l+papf*. (2.20) 

This standard boost takes the standard time-like .D-vector (1, 0, 0, . . .) into the unit normal- 
ized D- vector with spatial components given by p a , 

a( ^(s)X;) A ' (2 - 2i) 

A general vielbein can now be written as the standard boost of an upper triangular vielbein 

E A = Mv) A E B ( Nl + N%apa ^ + ^ W + ^^ ) (2 22) 

" V <Va e,\K + ^iP b P a ) J' 

This is simply a reparametrization of a general vielbein, one which will be particularly 
convenient for the Hamiltonian analysis. There need not be any gauge or Lorentz invariance 
to do this. The D 2 components of the general vielbein are now parametrized by the D 
components of N and N l , the (D — l) 2 components of the spatial vielbein e^, and the D — 1 
components p a . 

We now express the Einstein-Hilbert term in terms of this decomposition. The Einstein- 
Hilbert term is invariant under local Lorentz transformations. Therefore, when we plug in 
the vielbein as parametrized in (2.22), all the p a dependence drops out. Thus we can evaluate 
the Einstein-Hilbert action using the upper triangular ansatz (2.17). 



2 This is analogous to the standard boost used to define single particle states in Lorentz invariant quantum 
theory. See for instance chapter 2 of [33]. 
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The Hamiltonian formulation of GR in upper triangular vielbein form is reviewed in 
Appendix B. The result is that the Einstein-Hilbert kinetic term can be written in the form 
of a constrained Hamiltonian system on the 2d 2 dimensional phase space consisting of the 
spatial vielbein components e^, and their canonical momenta ir l a , 

Seh = J d D x (V a e 4 a -NC- N% - l -\ ah V ab \ . (2.23) 

Here C(e, ir), Ci(e, n) are the usual diffeomorphism constraints of GR, whose Lagrange multi- 
pliers are the lapse and shift. They depend only on the spatial vielbeins and their conjugate 
momenta. In addition, we have d(d — l)/2 primary constraints V a b{e, tt) = ej[ a 7r J fc ,, responsi- 
ble for generating the residual spatial local Lorentz rotations of the upper triangular vielbein, 
along with their Lagrange multipliers \ ah . 

For a theory with M spin-2 fields, each spatial vielbein e(i\ a gets canonical momenta 
7T(if a . Each of the Einstein-Hilbert terms is separately Lorentz invariant, and so will not 
depend on the P(/) a 's, 

Seh = f d D x ]T Uij a e {I)i a - N (1) C W ~ Nrf C {I)i - ^\ {I f V Wab \ . (2.24) 

Here C(i)(e,ir), C{i\{e, tt), are the diffeomorphism constraints of GR, one for each of the 
diffeomorphisms of the M Einstein-Hilbert terms, whose Lagrange multipliers are the lapses 
and shifts. In addition, there are the M sets of additional primary constraints, 

'PWab&n) = e Wi[a -K { i) b] , 

responsible for generating spatial local Lorentz rotations of each of the spatial vielbeins, 
along with their Lagrange multipliers A(/™ . 

Let's now consider adding a general interaction term (2.6). We assume the theory 
graph for our interaction term is connected. This is no loss of generality, since the various 
connected islands of a disconnected graph do not interact with each other and can be treated 
independently. The proposed potential terms thus have only one overall Lorentz invariance 
so we are not free to choose a gauge for each vielbein. For each vielbein we use the general 
parametrization (2.22), 

A I N mn + N w l e {I)l a m a N (mi f + N (I fe (I)i b (S b a + ^Pd^df) \ 
E<-% = I h • (2.25) 

\ Z(Di a W) a e ( /) 4 b (V + ^jPWbPdf) J 
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We may use the overall local Lorentz invariance to set one of the p a, s to zero, say paf = 0. 
This leaves a residual overall local spatial rotation invariance. 

In the Hamiltonian (2.24), the spatial vielbein is unconstrained, and there is a spatial 
local Lorentz rotation left over as a gauge symmetry, enforced by the primary constraints 
V(i) ab . Alternatively, we can if we wish eliminate the spatial local Lorentz symmetry com- 
pletely and work with a constrained spatial vielbein. To do so we decompose e^" as a 
rotation of some constrained spatial vielbein, 



em,- 



° = (e^>)W- (2-26) 



Here the qo) are anti-symmetric matrices 5 ac qw c = —S bc q(if c and (e 9(7) ) a fe is the usual matrix 
exponential. For each J, this is simply a general SO(d) rotation, parametrized by the 
d(d — l)/2 components of g(i> By e^ 6 , we denote a spatial vielbein satisfying, for each 
J, d(d — l)/2 constraints we specify shortly. Thus each constrained vielbein e^ contains 
d(d + l)/2 components. 

The constraints satisfied by the e^ 0, are chosen to be the following: 

■7^^e Wl d bc = -7^— c ^\ d ha • (2.27) 

where U is the full interaction potential (2.6). This symmetrization condition puts \d(d— 1) 
constraints on each eg^", eliminating ^d(d — 1) independent components, which have been 
transferred into the q(i). 

Using the parametrization (2.26), the Einstein-Hilbert kinetic terms will be indepen- 
dent of the q(i) b a , as they are independently invariant under spatial rotations, so the qa) b a 
are non-dynamical and will appear only in the potential term. The Einstein-Hilbert kinetic 
terms will now read 

S EH = / d D x Y, (W^° - N {1) C w (e, tt) - Ntf C (I)l (e, *)) . (2.28) 

This looks just like (2.24), except the spatial vielbeins are the constrained ones, and the 
canonically conjugate momenta 7t(if a contain only d(d+ l)/2 independent components since 
they are the conjugate momenta for the e^" alone and not the full spatial vielbein. There 
are no primary constraints V(i) ab necessary. 
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Consider now the qo) equations of motion, obtained from variation of the potential U 
with respect to the g(i), 

6U = ^ 5q ^ = ° • (2 ' 29) 

Since the matrix 5 ac q(i) b c is anti-symmetric, the equations of motion read 

9U S b ^^Ls ba . (2.30) 



dq { i) a b ' dq { iy b 
Rewriting this in terms of the decomposition (2.26) gives 



-^r [e _9CJ) ] b e Wl e (e^) d e 5 dc = — r " ' n " 



[e-**] \ e (1)l e (e^) d e 5 dc = -^ [e^«] \ e Wj e (e^) d e 5 da . (2.31) 



Consider expanding this in powers of qw. Due to the condition given in (2.27), we see that 
the zeroth order term vanishes, so the equations of motion are solved by 

qm\ = 0. (2.32) 

Thus the qo) are auxiliary variables, and can in this way be eliminated from the action by 
their own equations of motion. We may work only in terms of the constrained spatial vielbein 
td)"" by simply replacing t(i)°- — \ e^." in the potential. 

We note that there are constraints analogous to (2.27) for the full spacetime vielbein 
E A , which we consider in Appendix C. 

2.3 Ghost-freedom 

We now demonstrate the existence of M — 1 Hamiltonian constraints of the M- vielbein theory. 
This will guarantee the right number of degrees of freedom to describe one massless graviton 
and jV — 1 massive gravitons in D dimensions. 

To identify the primary constraints, let us consider now the interaction term (2.6), 
with the vielbeins given by the parametrization (2.25). We observe the crucial fact that 
both E(i) ° and E^ are linear in the lapses N(i) and shifts N(i), while E^ and E^ are 
independent of the lapses and shifts. Thus, due to the structure of the epsilon tensor, the 
interaction terms (2.6) are manifestly linear in all the lapses and shifts and can be written 

in the form 

AT 

U = J2 ( N ^ C ^ m + N ^ C Vk) ■ (2-33) 



7=1 
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where C(i) m (e,p), C(/™(e,p) are functions of the spatial vielbeins e^" and the boost vectors 
P(if. The "m" superscript denotes that these terms are coming from the mass (interaction) 
term. 

The full action is given by 
S = Jd D xf2 (W<M° - N {1) [C w + C (I) m ] - Ntf [C Wl + C«T } - ^A W °% )( J . (2.34) 

We now use M — 1 of the M shift equations of motion to solve for the M — 1 remaining p(i) a 
variables (recall that we have used the overall local Lorentz invariance to rotate away P(i) a ), 

Qi k (e,n)+C ( i?(e,p) = =► p w a = p {I) a (e,n) , I = 2,...,Af. (2.35) 

We have thus eliminated all the p(i) from the action in favor of the e^i) and 7T(i> After doing 
this, the action remains linear in the M lapse variables N(i) and the one remaining shift 
variable, Nof, 

S= [d D xlj2 (^a^\ a ~ Nu) [Cm + C w m ] - \\D ab V Wa )j - N<j [C(d 4 + C (1) ° 

'(2.36) 

There is a lapse and shift constraint with No.) and No) as multipliers, which will (after 
mixing with other constraints) generate the overall unbroken diffeomorphism symmetry of 
the theory. There are Af x \d(d — 1) spatial local Lorentz (i.e, rotation) constraints. One 
combination, the overall sum Yli^ii) &> w ^^ generate the unbroken overall spatial local 
Lorentz invariance. The other combinations will lead to secondary constraints. 

Finally, and most importantly, we see the presence of the additional TV— 1 primary con- 
straints with N(2), . . . , N(Af) as multipliers. These are the constraints responsible for removing 
the potentially Boulware-Deser ghost-like modes. As a theory with A^ interacting metrics 
will have a spectrum consisting of one massless graviton and M — 1 massive gravitons, there 
are M — 1 potential Boulware-Deser ghosts, one for each massive graviton. Here we see that 
for the interaction terms (2.6), there are just the right number of extra primary constraints 
to vanquish all of them. 

Let us count degrees of freedom. The phase space starts with J\f x 2d 2 degrees of 
freedom of the spatial vielbeins and their conjugate momenta. There is one overall spatial 
local Lorentz invariance generated by the first class constraints ^2];V(i) ab , thus removing 
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2 x d(d— 1)/2 dimensions of phase space. The other N '— 1 combinations of Lorentz generators 
will lead to secondary constraints, which together with the primary constraints will generate 
a second class set, removing an additional 2 x (A/" — l)d(d — l)/2 dimensions of phase space. 
We can see that these secondary constraints are indeed generated and that they form a 
second class set by using the parametrization (2.26). With this parametrization, the phase 
space contains only the constrained vielbein e(i) and its canonical momenta, for a total of 
2 x J\fd(d + l)/2 phase space dimensions. (See also Appendix B.) 

There is one overall diffeomorphism, the generators of which enforce first class con- 
straints, thus removing 2 x (d + 1) dimensions of phase space. Finally, we have the A/" — 1 
additional primary constraints found above. Though we will not show it here, we expect 
these constraints to each generate a secondary constraint, with which they should form a 
second class set, thus removing a further 2 x (A/"— 1) degrees of freedom. Adding everything 
up, we find that the physical phase space should have 

2 [\d{d - 1) - 1] + 2{N - 1) [\d{d + 1) - 1] , (2.37) 

dimensions. This corresponds exactly to one massless spin-2 field and A/" — 1 massive spin-2 
fields and their canonical momenta, with no extra ghosts. 

It should be emphasized that we have only shown the existence of the primary con- 
straints necessary for eliminating the ghosts. To complete the proof that these theories are 
ghost-free, the secondary constraints (which arise by demanding that the primary constraints 
be preserved in time) must be computed. To show that there are not too few degrees of free- 
dom, it must also be argued that there are no further constraints, and that the primary and 
secondary constraints together form a second class set. However, it is hard to imagine that 
the secondary constraints would be absent, as they have been shown to exist in the bi-metric 
case (in the metric formalism) in [11] and, what's more, their absence would indicate the 
existence of some kind of half degree of freedom which would not be consistent with Lorentz 
invariance in dimensions greater than 2. (Though unlikely, it is possible that with multiple 
metrics, half degrees of freedom could pair up into full degrees of freedom.) 

The arguments here rely only on the fact that the Hamiltonian written in terms of the 
ADM variables is linear in all the lapse and shift variables, and goes through for the general 
interaction term (2.6). This includes any kind of theory graph, including tree graphs, loop 
graphs and graphs including tri-vertices, tetra-vertices, and beyond. The argument is fully 
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non-linear, and is valid to all orders in the fields and beyond any decoupling limits, such as 
those considered in [2]. 

3 Bi- gravity 

In this section we show that, for the case of A^ = 2 dynamical spin-2 fields, the vielbein 
theories introduced above are dynamically equivalent to the metric ghost-free bi-gravity 
theories, studied in [18]. The metric bi-gravity theory is obtained by promoting the flat 
reference metric rj^ of dRGT massive gravity to a general metric / M „ and allowing it to be 
dynamical by adding an Einstein- Hilbert term for f^, with its own Planck mass Mf, 



S = d u x 



M°~ 2 __ Mf' 2 ,— m 2 M?- 2 _^ ' 



n=0 



(3.1) 



where Mf~ 2 = {1/M°~ 2 + l/Mf~ 2 ) 



There are D + l different symmetric polynomials, and hence D + l different parameters 
in the interaction term. The first symmetric polynomial 1/— g Sq{^ g~ l f) = y/—g and the D- 
th polynomial \J—gSD(\/g~ l f) = y/—f are cosmological constants for g and /, respectively. 
Each polynomial gives a tadpole when expanded around flat space: g^ v = rj^ + -^-g^ v and 
ffM> — Vnu + M~fw If we demand that flat space is a solution for both metrics (i.e., no 
tadpoles), then we must take 

f^ Q k\{D-k)\ ' f-^ (k -1)\{D -k)\ K } 

For D = 4 this means 

A> = - (3/?i + 3/3 2 + p 3 ) , A = - (A + 3/3 2 + 3/3 3 ) . (3.3) 

With this choice, expanding to quadratic order gives the Fierz-Pauli term [34] for the fluctu- 
ation jj-g^v — ~M~.f\iv The orthogonal fluctuation -^-g^ v + jfffiu is absent from the potential. 
Thus this theory propagates precisely one massive spin-2 field and one massless spin-2 field 
around flat space. We can absorb one further coefficient by taking m to be the mass of the 
massive spin-2 and setting 
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For D = 4, this gives 

f3i + 2/3 2 + /? 3 = 8 . (3.5) 

In -D dimensions, the bi-gravity theory will have D — 2 free parameters after tadpoles are 
eliminated and one coefficient is absorbed into the mass. 

Despite the asymmetric appearance, the potential in (3.1) does not actually favor one 
metric over the other, since we have the property 



-g 



SniVa^f) = y^fSn^i^f^g) . (3.6) 



If we wish, we can impose a Z 2 symmetry under the interchange g^ v -H- f^ u , by setting 
M g = Mf and (3 n = (3 D . n . 

The two Einstein-Hilbert terms are each invariant under a separate diffeomorphism 
symmetry, 

9^(x) -> dxii dxv g a p (/ ( i)(x)) , /^(a;) -)• ^7^7^ (/«(*)) ■ ( 3 - 7 ) 

The mass term breaks this down to the subgroup of diagonal diffeomorphisms with f{xf = 
fof = / M , which acts the same way on both metrics, 

df a dfP df a df 13 

9 ^ X) ~* d^d^ 9a/B (/(x)) ' f ^ x) ""* d^d^ 9a0 (/(x)) ' (3 ' 8) 

Naive counting arguments would suggest that a generic theory of two interacting spin-2 fields 
with only one diffeomorphism invariance will propagate 8 degrees of freedom for D = 4: one 
massless spin-2, one massive spin-2 and one scalar ghost. However, it was shown in [11, 18] 
that the theories (3.1) propagate only 7 degrees of freedom at the full non-linear level, 
consistent with one massless and one massive spin-2 alone. 



3.1 Vielbein formulation of bi-gravity 

We show now that the metric bi-gravity theory has an equivalent vielbein formulation. We 
introduce two vielbeins, one for each metric, 

9\u> = E(i) E mi/ r] AB , f^ u = Em E( 2 ) u r] AB . (3.9) 
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When the substitution (3.9) is made into the Einstein-Hilbert part of the action, it becomes 
invariant under two separate local Lorentz transformations, 

E Wn ~> A <!) B E Wfi ' E <-% -* A ( 2 ) B E <-% > ( 3 - 10 ) 

in addition to the two diffeomorphisms, under which the vielbeins transform as one-forms, 
Ev£(x) -> |^W (Mx)) , E (2)f f(x) -> ^W (/raO*)) . (3.11) 

We will show that the following action, which consists of the Einstein-Hilbert terms in 
vielbein form and an interaction term written using wedge products of the vielbein one-forms 
E A = E A dx^, is dynamically equivalent to (3.1) in terms of the variables (3.9), 

S = — |— / d D x (detE m )R[E { i)} + —$— / d D x (det E (2) ) R[E (2) } (3.12) 

m 2 M D ~ 2 f D R 

£- / T uj : n ~e AlA2 -A D ^ M A • • • A E (2 A A EWW A • • • A E W A ° . 

8 / *-^ n\[D — n)\ 

J n=0 V ' 

The mass term breaks the two diffeomorphism and local Lorentz symmetries down to the 
diagonal subgroup where A(x) B = A( 2 ) B and /(i) M = /(2) M - 

Using the relation (2.12), the potential can be written in terms of the symmetric 
polynomials (defined in Appendix A.l), 

D 

J2 Pn (det E m ) S n (E m - l E m ) . (3.13) 

n=0 

Our first step is to show that the action with this potential, is in fact equivalent to the same 
action with the additional constraint that the following product of vielbeins and vielbein 
inverses be symmetric with respect to the Minkowski metric 3 . In matrix notation, 

E^E^r] = 7] (E(i ) - 1 Ev)) T . (3.14) 

To see this equivalence, parametrize one of the vielbeins, say the first vielbein E(i), as a 
generic Lorentz transformation times a constrained vielbein E(i), chosen to satisfy (3.14), 

E w = E w e u , E^ 1 E m = r) (e^ 1 £ (2) ) • (3.15) 



3 Note that, taking the inverse of both sides of (3.14), we obtain E&) E^r] = r)(E&) 1 E(i)) , so this 



condition is in fact symmetric under lf>2. 
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We have written the Lorentz transformation as the exponential of a matrix u which is anti- 
symmetric with respect to i], 

r/uj = — u T r] . (3.16) 

Equation (3.15) is nothing but a parametrization of the general vielbein - we have packaged 
the D 2 components of the general vielbein into the D(D + 1)/2 components of a constrained 
vielbein and D(D — l)/2 components of a Lorentz transformation. 

Now, we will see that the D(D — l)/2 variables in u are auxiliary variables, and that 
their equations of motion set u = 0. Plugging the decomposition (3.15) into the action (3.12), 
the parameters in the Lorentz transformation only appear through the potential, since the 
Einstein-Hilbert term is separately Lorentz invariant in each vielbein. The potential takes 
the form 

D 

J2 /Wet E w ) S n (e-"Em X E { ^ . (3.17) 

n=0 

Now vary with respect to u. Consider first the lowest order terms in u which contain no 
powers of u beyond 5u. Due to the form of the symmetric polynomials, the only terms 
that appear at lowest order are traces of 5u with powers of the matrix E(t) Ep). Since 
E(i) E(2) is symmetric and 5u antisymmetric (both with respect to rj), and since 5u appears 
only linearly, the lowest order expressions in u vanish. Thus the equations of motion start 
linearly in u, and are solved by 4 

co = 0. (3.18) 

Thus the action with unconstrained vielbeins is dynamically equivalent to the action with 
constrained vielbeins. We may plug u> = into the potential (3.17), 

D 

J2 /3n(det E w )S n (E w ~ l E (2) ) . (3.19) 

n=0 

Let us now relate the vielbein potential to the metric potential. In matrix notation the 
metrics are given by g = E^r] E(i) T and / = Ew rj E^f . Therefore, 

<T7 = (Em-y r]- 1 E w - 1 EvwEp? . (3.20) 



4 Thcre is also the possibility of having non-trivial solutions of the ui equations, in which case there would 
be more than one branch of the theory. In this case, the uj = branch of the vielbein theory would be 
equivalent to the metric theory, and the other branches may not be. 
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Using the parametrization (3.15), along with the symmetry property of E{\) gives g 1 f = 
E(i) E{2) E(i) E(2) , or 

V / F T 7= (E^E^Y . (3.21) 

Using the properties (A. 8) and (A. 9), we may write 

(det £ (1) ) S n (E w ' l E w ) = ^/-detg S n (^g^]) . (3.22) 

We see that the vielbein bi-gravity theory is dynamically equivalent to the metric bi-gravity 
theory, as claimed. 

Note that we have only proved the equivalence of the metric and vielbein formulation 
for theories with one bi-gravity interaction. When constructing multi-gravity theory graphs 
out of the bi- vertex as in figure 2, the equivalence of the vielbein formulation and the metric 
formulation does not appear to always hold. In Appendix C we argue that the equivalence 
holds as long as the theory graph is a tree graph, i.e., contains no closed loops. If the theory 
graph contains loops however, the equivalence appears to break down. This is directly 
related to the form of the constraints used to pass from the vielbein formulation to the 
metric formulation (3.14). Essentially, the loop graphs alter the form of these constraints. 

In the vielbein formulation, the proof of ghost-freedom is trivial to extend to all theory 
graphs. Thus we expect graphs such as figure 2 to be ghost-free when the interactions are 
given in the vielbein formulation. For the metric theory however, it is not known whether 
the proof of ghost-freedom extends to the case of theory graphs with closed loops. The 
multi-gravity theory represented in figure 2, or, for example, a "triangle" theory such as that 
studied in [35], potentially contains ghosts when formulated in terms of metrics rather than 
vielbeins. (See Appendix C for more on this point.) 

4 dRGT massive gravity 

We now turn to the case of massive gravity. We treat this case separately not only because of 
specific interest in massive gravity but also because these theories do not inherently possess 
the overall diffeomorphism invariance of the multi-metric theories. Thus in this section 
we will repeat the Hamiltonian analysis specifically for massive gravity. We will see that 
demonstrating the existence of the primary constraint which eliminates the Boulware-Deser 
ghost is much easier in vielbein variables. 
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In D spacetime dimensions the action for ghost-free dRGT massive gravity is 

,2 D 



s=*£lu>, 



-g 



tf-T-E'WVS^) 



n=0 



(4.1; 



The massive gravity theory depends on the dynamical metric g^ and the fixed background 
metric rj^. The Einstein-Hilbert part of the action (4.1) is invariant under diffeomorphisms 

df a df 



9nu[x) ->■ 



;9ap (f(x)) 



(4.2) 



dx^ dx l 

but the mass term breaks this symmetry, due to the appearance of the background metric 

As in the bi-metric case, there are D + 1 different symmetric polynomials, and hence 
D + 1 different parameters in the mass term. For massive gravity, the D-th symmetric 
polynomial is just a/— det r] = 1, so it doesn't contribute to the equations of motion. Thus 
there are only D free parameters. We can ensure that flat space is a valid solution by 
demanding 



D 



«E 



fc=0 



k\(D-k)\ 



(D-iy-Yl 



' (k-l)\(D-k)\ ' 



(4.3) 



For D — 4, this gives 

0o = - (3/?i + 3/3 2 + h) . (4.4) 

Then, expanding to quadratic order gives the Fierz-Pauli term [34] for the fluctuation h^ v = 
9inv — Vt+v Thus this theory propagates precisely one massive spin-2 field around flat space. 
We again absorb one further coefficient by taking m to be the mass of the massive spin-2 
and setting 

_£>! V" £* + (D-2V V ^ 



Q k\(D-k)\ v '■ ^ (k - 2)\(D - k)\ 



(4.5) 



For D = 4, this gives 



/3 1 + 2/3 2 + /3 3 = 8. (4.6) 

In D dimensions, the theory has D — 2 free parameters in addition to the mass. 
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4.1 Vielbein formulation of massive gravity 

Our goal now is to show that dRGT massive gravity in D dimensions (4.1) has a dynamically 
equivalent vielbein formulation 5 . We introduce the vielbein fields E A (x), 

g^ = E A E u B VAB , (4.7) 

along with the following action, which consists of the Einstein-Hilbert term in vielbein form, 
and potential terms written using wedge products of the vielbein one-forms E" 4 = E A dx^, 
and unit one-forms which can be thought of as vielbeins for the flat background metric 
1 A = 5 A dx^ 

m d ~ 2 ( r 

S = — p — \ d D x det{E)R[E] 

~ X / ^ n\(D- n)\ ~ eAlA *-" AolAl A --- AlAnA EAn+1 A • • ' A eA °) ■ ( 48 ) 

By pulling out a determinant as in (2.12), the mass term can also be written in terms of the 
symmetric polynomials of the matrix E A , 

1 



n\{D-n)\ 



€ AiA2 ... Ad 1 Ai A . . . A l An A E An+1 A . . . A E Ad = {detEd D x)S n {E- 1 ) . (4.9) 



The action (4.8) is an action for D 2 variables, whereas (4.1) is an action for D(D + 1)/2 
variables. Furthermore, (4.8) has no gauge symmetry, since both diffeomorphisms and LLI 
are broken by the mass term. Nevertheless, we will show that (4.8) is dynamically equivalent 
to (4.1). The logic proceeds much as it did for the bi-gravity case. 

We show that the action (4.8) is dynamically equivalent to the same action with the 
additional constraint that the vielbein be symmetric with respect to the Minkowski metric, 

Er] = r]E T . (4.10) 

We parametrize the general vielbein as a constrained vielbein E satisfying (4.10), times a 
generic Lorentz transformation, 

E = Ee". (4.11) 



3 Vielbein massive gravity is also considered in [21]. 
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Again, the Lorentz transformation is written as the exponential of a matrix u which is 
anti-symmetric with respect to 77, 

T]LO = -u T Tj. (4.12) 

The w's will not appear through the Einstein-Hilbert term, since it is invariant under local 
Lorentz transformations. 

As in bi-gravity, the D(D — l)/2 variables in u appear only through the mass term and 
are auxiliary variables. Their equations of motion will set u = 0. Plugging the decomposition 
(4.11) into the action (4.8), the parameters in the Lorentz transformation appear through 
the mass term, 

D 

^/3„(det E)S n (e-»E- 1 ) . (4.13) 

n=0 

Now consider the u equations of motion. Vary with respect to u and consider the equations 
of motion in powers of u. The lowest order terms contain no powers of u beyond the variation 
5u. The only terms that appear at lowest order in the symmetric polynomials are traces of 
5u with powers of E~ 1 , Since E' 1 is symmetric and 5u antisymmetric (both with respect 
to rj), and since 5u appears only linearly, the terms in the equations of motion linear in u 
vanish. Thus the equations of motion start linearly in u, and are solved by 6 

u = 0. (4.14) 

We may plug the solution (4.14) into the potential (4.13), giving 

D 

J2(3n(detE)S n (E- 1 ). (4.15) 

n=0 

Thus the action with unconstrained vielbeins is dynamically equivalent to the action with 
vielbeins constrained to satisfy (4.10). 

To relate the potential (4.15) to the metric potential we write g = E rj E T so that 

g- 1 r ] = (E- l ) T r]- 1 E- 1 <q . (4.16) 

Using the parametrization (4.11) and the symmetry property of E~ x , we see that 

y^={E- 1 f. (4.17) 



6 As in the bi-metric case, there is the possibility of additional branches if there are non-trivial solutions 
of the to equations. If these exist, then the equivalence between vielbein and metric theories may only hold 
for the trivial branch. 
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Thus we can write, using the property (A. 8) of the symmetric polynomials, 



(det£) S^E" 1 ) = ^-detg S n {^f^q) . (4.18) 

We see that the vielbein massive gravity is equivalent to dRGT massive gravity. 

4.2 Ghost-freedom 

We've seen that we can write the action for dRGT massive gravity using the D 2 components 
of an unconstrained vielbein as variables. Now, by choosing a different parametrization for 
the vielbein, we will see that it is almost trivial to identify the primary constraint which 
eliminates the Boulware-Deser ghost. 

We again perform a d + 1 decomposition of the general vielbein as in (2.22), an upper 
triangular vielbein rotated by a standard Lorentz boost parametrized by p a , 

e a = ( N<y + Wefpa Np a + N\ b (6 b a + ^p b p a ) \ 

" ' V e i a p° e M° + ^iPbP a ) J ' 

The p a do not enter the Einstein-Hilbert term, since it is Lorentz invariant. 

The mass term (4.9) is not invariant under local Lorentz transformations, so there will 
be explicit dependence on the p al s. Note, however, that the lapse and shift, N and N l , only 
appear in the components E ° and E b , and that they both appear linearly. Due to the 
epsilon tensor in the mass term, there will never be more than one component E ° or E b in 
any term, and so the entire interaction term is manifestly linear in both the lapse and the 
shift iV and N l . Thus we can write the mass term in the form 

U = NC m (e, p) + NX™{e, p) + H(e, p) . (4.20) 

As the Einstein-Hilbert Hamiltonian is also linear in the lapse and shift, these remain La- 
grange multipliers in the full massive theory, enforcing the constraints 

C(e,7r)+C m (e,p)=0, ^(e, vr) + Cf(e,p) = . (4.21) 

Note that in the metric formulation, the lapse and shift do not automatically appear in this 

way. 
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The p a, s are auxiliary variables, since they appear appear only through the mass term, 
with no derivatives. We can eliminate them by solving the N* constraints, 

&(e, tt) + C(e, p) = =► p a = p a (e, vr) . (4.22) 

The action now takes the form 

S = J d D x (n\e l a -n(e,p(e,n))-^\ ab V ab (e,n)-N[C(e,n)+C m (e,p(e,7r))]\ . (4.23) 

We see explicitly the presence of the extra primary constraint, enforced by the shift N, 
responsible for removing the Boulware-Deser ghost. As shown in [11], there will be an addi- 
tional secondary constraint which arises by demanding that the primary constraint enforced 
by N be conserved in time. 

Our theory has no gauge symmetry - no diffeomorphism symmetry or local Lorentz 
invariance - so there will also be secondary constraints which come from demanding that 
the primary Lorentz constraints V a b be preserved in time. These combine with the primary 
constraints into a second class set. Alternatively, we may fix the local rotation invariance by 
parametrizing the spatial vielbein as the rotation of some constrained spatial vielbein (2.26) 
as was done in the case of multi-metric gravity in subsection 2.2 as well as in Appendix B. 

Let us now count degrees of freedom: the phase space has 2d 2 variables, the d 2 spatial 
components of the vielbein e^, and their d 2 canonical momenta 7r* a . There are then second 
class constraints restricting the phase space: 2 x ^ ~ - constraints coming from the local 
Lorentz constraints and their secondary constraints, and 2x1 special constraints coming 
from the lapse multiplier and the corresponding secondary constraint. This leaves a 

2 (*^> - l) (4.24) 

dimensional physical phase space, just right to describe the degrees of freedom of a traceless 
symmetric tensor and its canonical momenta, i.e. a massive spin-2 graviton, with no extra 
Boulware-Deser modes. 

Given a background solution to the bi-gravity theory (3.1) where one of the metrics is 
Minkowski, we can recover the massive gravity theory (4.1) by sending the Planck mass asso- 
ciated with the Minkowski metric to infinity, after canonically normalizing the fluctuations. 
This can be done for any background other than Minkowski, and the result will be ghost-free 
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massive gravity propagating on that background. In fact, many different theories containing 
only massive gravitons, around various and even multiple backgrounds, can be obtained by 
taking the appropriate M{i) — > oo limits of the general multi-metric theory. By arguments 
similar to those in this section, it can be seen that all these theories are ghost-free. 

5 Discussion 

In this paper, we have reformulated the recently uncovered theories of massive gravity and 
bi-gravity in terms of vielbein variables. We find that with this choice of variables, the 
theories become much more natural and transparent. In place of the unwieldy polynomials 
and matrix square roots of the metric formulation, the vielbein formulation simply contains 
wedge products of the possible combinations of the vielbeins. In vielbein form, the additional 
primary constraints which signal the absence of the Boulware-Deser ghost are manifest, to 
all orders in the fields and beyond any decoupling limits. 

The vielbein formulation has allowed us to extend the ghost-free interactions to in- 
teractions between multiple metrics. The natural extension is to allow for arbitrary wedge 
products between all possible vielbeins, leading to the general interaction (2.6). We have 
seen that the proof of the existence of additional primary constraints extends easily to this 
general interaction. 

We conjecture that the interaction (2.6) is in fact the most general ghost-free zero- 
derivative interaction among multiple vielbeins, and thus the most general such interaction 
among massive and massless spin-2 fields. This interaction can be though of as a kind of 
multi-metric generalized cosmo logical constant, or a lowest-order generalized multi-Lovelock 
invariant' [37, 38]. We will expound on this relation in future work. 

This construction opens up a wide parameter space in which to search for phenomeno- 
logically interesting or promising models, and in which to search for theories which avoid 
potential problems, such as superluminality around non-trivial solutions, and strong coupling 
around self-accelerating solutions. It has been argued that for a single massive graviton of 
mass m in D = 4, the highest unitarity bound which is possible is the rather low cutoff 



7 Note that the addition of the usual single-metric Lovelock invariants to the multi-gravity theories does 
not disrupt the proof of ghost-freedom [36] . 



A3 = (Mpm 2 ) [31]. It is possible that within these interacting multi-metric theories, there 
are examples with higher cutoffs. 

Finally, we comment on the possible coupling of multiple metrics to matter. If matter 
is minimally coupled to only a single metric, then clocks and rulers will measure distances as 
determined by that single metric. This kind of coupling will not re-introduce the Boulware- 
Deser ghost, as it maintains the same symmetries of the Einstein-Hilbert term. However, 
more general ghost-free couplings might exist that involve more than one metric and maintain 
an overall diffeomorphism invariance (we might call this multi- minimal coupling) , though this 
can naively be expected to violate the equivalence principle. 
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A Symmetric polynomials 

In this appendix, we define the various matrix polynomials used throughout the paper, and 
describe some of their properties. 

A.l Elementary symmetric polynomials 

Given a D x D matrix M A B , we define the elementary symmetric polynomials, for < n < D, 
Sn(M) = n[ ^_ n)[ e Al A 2 ...A D l B ^- B ° M% ■ ■ ■ M\8 A £l ■ ■ ■ 8% , (A.l) 
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or equivalently 



M% ■ ■ ■ M a S J A bZ ■ ■ ■ S A £ d Zmm...a d = nKD D[ ^ S n (M)~e BlB2 ... BD . (A.2) 

In terms of traces of products of M, the first few are 

S (M) = 1, 

S^M) = [M], 

S 2 (M) = ^([M] 2 -[M 2 ]), 

S 3 (M) = i([M] 3 -3[M][M 2 ] + 2[M 3 ]) , 

s 4 (m) = ^([m] 4 -6[m] 2 [m 2 ]+8[m][m 3 ]+3[m 2 ] 2 -6[m 4 ]) , 

; (a.3) 

The D-th symmetric polynomial is the determinant, 

S D (M)=detM, (A.4) 

and the higher symmetric polynomials are defined to vanish identically, 

S n (M) = for n> D. (A.5) 

If M is diagonalizable, the symmetric polynomials are the symmetric polynomials in 
the eigenvalues. If we label the eigenvalues (including degeneracy) \a, A — 1,---,D, then 



S (M) -- 


= 1, 


S 1 (M) -- 


= 5>. 

A 


S 2 (M) = 


= 2_^ ^A^B , 




A<B 


S 3 (M) -- 


/ _, A^AsA. 



c 

A<B<C 



S D {M) = AiA 2 ---A D . (A.6) 



30 



The symmetric polynomials may be obtained from expanding out the following deter- 
minant, in powers of e, 

D 

det (1 + eM) = J2 £ n Sn{M). (A.7) 

?1=0 

Transposing a matrix leaves its symmetric polynomial unchanged, 

S n (M T ) = S n (M), (A.8) 

and given the symmetric polynomial of a product of two D x D matrices M and N, we may 
cyclically permute the argument, 

S n (MN) = S n (NM). (A.9) 

A. 2 Generalized symmetric polynomials 

Define the generalized symmetric polynomials for multiple DxD matrices M^\ M^ 2 \ • • • , M( D \ 
as follows, 

S (M« M< 2 \ • • • , MM) = l AlM ... AD &*•"*» mV\mV>\ ■ ■ ■ M^X , (A.10) 

or equivalent ly 

M (1) > (2) t • • • M^%l AlM ... Ao = ^S (M«, M< 2 \ -.., M<">) ~e BlB2 ... BD . (A.11) 

We can calculate 5 (M' 1 ', M^ 2 \ • • • ,M^ D )) by writing all D\ of the top to bottom con- 
tractions of the indices of the M's, with a sign for straightening out the contraction. For 
example: 

• for D = 2, there are two possible contractions, and we have 

S(X,Y) = [X][Y]-[XY], (A. 12) 

• for D = 3 there are six possible contractions, and we have 

S(X,Y,Z) = [X][Y][Z] - [X][YZ] - [Y][ZX] - [Z][XY] + [XYZ] + [XZY], (A. 13) 
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• for D = 4 there are 24 possible contractions 

S(X,Y,Z,W) = [X][Y][Z][W] - [XYZW] - [XYWZ] - [XZYW] 

-[XZWY] - [XWYZ] - [XWZY] + [X] ([YZW] + [YWZ]\ + [Y] ([XZW] + [XWZ] 

+ [Z]([XYW] + [XWY]) + [W]([XYZ] + [XZY]\ - [X][Y][ZW] - [X][Z][YW] 
-[X][W][YZ] - [Y][Z][XW] - [Y][W][XZ] - [Z][W][XY] 

+[xy][zw] + [xz][yw] + [i^][yz] . 

(A.14) 

The elementary symmetric polynomials defined in (A.l) are the special case, 

S n (M) = ^^(M^, V^l). (A.15) 

n times D— n times 

B Hamiltonian formulation of GR in vielbein variables 

Here we review the Hamiltonian formulation of GR, in terms of upper triangular vielbein 
variables (2.17) (see also section 2.3 of [39]), 

*-(:?)-*-U'.)- 

The Lagrangian form of the action is 

M D-2 p M D-2 p 

S = — |— / d D x (det E) R{E) = — p — / d D x det(e) N (R[e] - K 2 + K i3 K %3 ) , (B.2) 

where surface terms are omitted and the extrinsic curvature is written in terms of vielbeins, 

Kij = ^ {e t a e ja + e/e m - V.iV, - V^) . (B.3) 

The covariant derivatives are with respect to the spatial metric gij, expressed in terms of the 
spatial vielbeins through g^ = efe^Sab- 

Now, we Legendre transform with respect to the spatial components of the vielbein, 
e°". The canonical momenta are 

v\ = ^ = M£" 2 det(e)e ja (K» - Kg*) . (B.4) 

32 



Multiplying this by an e^ and realizing that K^ is symmetric, we arrive at a set of primary 
constraints, 

V ah = e l[a 7i l b] = . (B.5) 

The V a b are anti-symmetric and represent \d(d — 1) constraints. 

We can invert for K^ in terms of ir t a , 

Ki 3 = —fr4 TT (V< - -TV—^A^k^l) ■ (B.6) 

M£- 2 det(e) \ D-2 bk l ) y ' 

From the definition of K^ , we have 

e {i a e j)a = NK ij +V {i N ]) . (B.7) 

Now we can calculate the Hamiltonian density, 

H = *\e t a - £ (B.8) 

M D ~ 2 
= — ^— det(e) [2e ja e, a (K^ - Kg**) - R[e] + K 2 - K^K i3 ] . (B.9) 

Using (B.7) and the symmetry of K^ to express the first term in terms of K's, 

2e ja e t a [ICJ - Kgv) = (K ij - Kh ij ) (2NK ij + 2V i iV i ) , (B.10) 



we find, after integrating by parts the terms containing N, 



t- 



M D ~ 2 
H = —?— [N det(e) (-R[e\ - K 2 + K ij K io ) - 2N t det(e) Vj (K ij - Kg ij )] . (B.ll) 

Now, using the expression (B.6), we have expressed the Hamiltonian in terms of the canonical 
momenta and vielbein. 

Adding Lagrange multipliers X ab for the primary constraints (B.5), the GR action is 
now in the form of a constrained Hamiltonian system, 

S= j d D x U^eS -NC- N t C l - l -\ ab V a A , (B.12) 

where 

iyfD-2 M D ~ 2 

C = — p — det(e) (R[e\ + K 2 - K* K %j ) , C = — ^— 2 det(e) V, (K** - Kg**) , V ab = e,^ 

(B.13) 
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Here's how the counting of constraints works in pure GR: the phase space is 2d 2 dimensional, 
containing the d 2 components of the spatial vielbein and their canonical momenta. The 
lapse and shift appear as Lagrange multipliers, enforcing the d constraints C l = and 
the 1 constraint C = 0. These are the diffeomorphism constraints associated to spatial 
reparametrizations and time reparametrizations, respectively. On top of that, we have the 
d(d — l)/2 additional primary constraints V a b, which generate the spatial local Lorentz 
symmetries of the upper triangular vielbein. These constraints are all first class, so each 

removes two dimensions from the phase space. We are left with a 2 ( ^—^ 1 ) dimensional 

phase space, just right for describing the degrees of freedom of a transverse symmetric 
traceless tensor mode, i.e. a massless graviton. 

Alternatively, we can directly fix the spatial local Lorentz invariance by parametrizing 
the spatial vielbein e t a as the rotation of some constrained spatial vielbein, 

e l a = (e«)W • (B.H) 

The q are anti-symmetric, 5 ac q b c = —S bc q a c and {e q ) a b is the usual matrix exponential, giving 
a general SO(d) rotation parametrized by the d(d — l)/2 components of q. By i b , we denote 
a spatial vielbein satisfying any d(d—l)/2 constraints which completely fix the local rotation 
invariance, so that if contains only d(d + l)/2 independent components. 

Inserting (B.14) into the action (B.2), the g's drop out, since the action is local Lorentz 
invariant. We may then Legendre transform with respect to if, obtaining canonical momenta 
Tx l a which contain only d(d + l)/2 independent components. The hatted variables will now 
identically satisfy (B.5), so there is no longer any primary constraint, and the Hamiltonian 
action takes the form 



S= d u x [ Ti\if -NC- NiC 1 , (B.15) 



with the constraints written as functions of the constrained hatted variables if, 7r* a . In 
this formulation, we do not need to worry about the Lorentz constraints of their secondary 
constraints, but we have to work with constrained spatial vielbeins. 



34 



C Lorentz constraints and vielbein/metric equivalence 



A theory of M interacting spin-2 fields should maintain the same number of degrees of free- 
dom when written in terms of vielbeins or in terms of metrics. A metric contains \D{D + 1) 
components while a vielbein contain D 2 components, giving an extra \D(D — \) components. 
The interaction terms containing M vielbeins considered in this paper have only one local 
Lorentz invariance, i.e. the overall invariance which rotates all the vielbeins together. This 
only accounts for the removal of the \D(D — 1) of one of the vielbeins. The remaining J\f — 1 
vielbeins must have their extra components removed in a different manner, and we will see 
in this Appendix how this comes about at the level of the Lagrangian. 

The extra components of the vielbein are removed by on-shell Lorentz constraints. To 
see this, let us consider first the usual Einstein-Hilbert kinetic term written in terms of the 
vielbein, 

S E h = \ I d D x (det E) R[E] . (C.l) 

This term is invariant under the local Lorentz transformation 

(C.2) 



(C.3) 



(C.4) 

Since the kinetic term is invariant under the local Lorentz transformation, the variation of 
the action under this transformation is zero, 

when 5E A = u A B E B . 







E* -> Eff : 


= A V/ 


Denote the infinitesimal Lorentz transformation by 






^B * Si 


+ uj a b , 


where co A B 


satisfies 










V CB co A B = - 


w B . 



Since u is an arbitrary anti-symmetric function of the spacetime coordinates, this 
implies that the following is identically zero off-shell 
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Consider now the multi-vielbein theory, E(i),E(2),. . ., with an Einstein-Hilbert kinetic 
term for each vielbein, as well as a potential term U(E(i),E(2), . . .) that mixes the E(i) but 
leaves precisely one overall Lorentz invariance intact. The E{i) equations of motion for this 

theory are 

1 5(det E (1) R[E (1) }) 5U(E W , E (2) , • • •) = Q , Q ?) 

2 SE(i )lM SEwp 

Multiplying both terms by E(i) t]cb and anti-symmetrizing, we find that the Einstein- 
Hilbert part vanishes due to (C.6), and the equations of motion imply 



5U(E W , Ew, . . ., Ir r 

A 



8E ( 



E Wu c r] CB -(A^B) = 0. {C.i 



n„ 



For each vielbein E(i), we therefore have one on-shell constraint, saying that the derivative 
of the potential times the vielbein is symmetric, 



5U(E(i), E(2), . . .) c SU(E(i), E( 2 ), . . .) r 

SE m * Em » » CB = -^K,„ 



-E(i) T] C B = r-— 5 E(i) 7] C A ■ (C.9) 



These constraints eliminate \D(D — 1) components of each vielbein. These constraints are 
invariant under the overall local Lorentz invariance. 

Note, however, that the sum of the constraints is identically satisfied, due to the one 
overall Lorentz invariance of the potential 

Thus for a (connected) theory with M vielbeins, there are in fact only M — 1 independent 
Lorentz constraints. The overall Lorentz invariance removes D(D — l)/2 components from 
the final vielbein. This means M x \D(D — 1) components of all the H vielbeins can be 
eliminated, leaving the same number of components as the ./V-metric theory. 

These constraints play a crucial role in relating the vielbein theories to equivalent 
metric theories. Consider the bi- vertex theory, given by (3.12). The potential is given by 
U = Ylm fin U n , for generic coefficients j3 n , where in matrix notation, the possible interaction 
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terms are 

U = det£ (1) , (C.ll) 

U x = det E m [E^E^ , 

U 2 = IdetEwdEw^E^-lEw-'E^Ew-'E^]) , 

U3 = g det E(i) ([-B(i) E(2)\ — 3[£<i) E(2)\[E(i) E(2)E(i) E(2)\ + 1\E(i) E{2)E{i) E{2)E(i) E(2)\] , 



There are two vielbeins, so we expect one independent Lorentz constraint. Using (C.9), 
it is straightforward to determine this constraint to be 

Ea^E^T) = (Eaf 1 E m ri) T . (C.12) 

(By taking the inverse of both sides, we see that this is equivalent to E(2)~ 1 E(i)T] = (E(2~ Ewq) 
so the constraint is in fact symmetric 1 -H- 2.) Note that it is precisely this constraint that 
was used in eq (3.14) in passing from the vielbein formulation to the metric formulation of 
bi-gravity. 

This constraint happens to be independent of the coefficients /3 n that appear in front 
of the various potential terms U n . This is a very nice property of the bi- vertex theory. 
As long as one considers only tree graphs of spin-2's that interact through the bi-vertex 
interactions given above, the constraints are of this simple form. Each line in the graph 
corresponds to one Lorentz constraint: in a tree graph with M vielbeins there are M — 1 
lines and M — 1 constraints. For a line that connects vielbein E(i) with vielbein E(j) one has 
the corresponding constraint E(i)~ 1 E(j)7] = (E(i)~ 1 E(j)T)) T '. This means that for every tree 
graph of bi-vertex interactions, there is an equivalent metric formulation of the theory, given 
simply by replacing E(if l E(j) with \/ g{i)~ l g(.j). 

As soon as one closes a loop in the graph, however, things get more complicated. The 
number of lines in the graph now exceeds the number of Lorentz constraints. To see what 
happens, consider the theory whose graph is a triangle, and where each side only contains 
the interaction U\ with a generic coefficient. The potential is 

U = an det E w [E w ~ l E^ + a 2 det E m [E^E^ + a 3 det E (s) [E^Em] , (C. 13) 

for generic coefficients a 2 , ct2, «3- The graph has three lines but only two Lorentz constraints, 
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given by enforcing the symmetry of the following two matrices, 

«i(det Ed)) £7 ( if x E { 2)T] - a 2 (det £7 (2 )) £7 (2 f x E^tj , (C.14) 

a 2 (det £'(2)) £7(2) -1 E^r) — a 3 (det £7( 3 )) E( 3 f E^r/ . (C.15) 

The constraints now depend on the details and coefficients of the potential. What's more, 
due to the form of the constraints, it's evident that the vielbein theory is no longer equivalent 
to the metric theory with the replacement E(i-r x E(j) — > \J gvr 1 gw- Thus, we suspect (but 
have not proven) that the theory of [35] has Boulware-Deser ghosts at the non-linear level. 

One can perform a similar analysis of the tri-vertex, 

U = det ^([^af^paf 1 ^)] - [E^E^E^E^]) . (C.16) 

There are again three vielbeins and thus two independent Lorentz constraints, given by 
enforcing the symmetry of the following two matrices, 

[£<i) .£7(2)] £7(i) £7(3) r\ — £7(i) £7(2)£7(i) £7(3) r\ , (C.17) 

[Ew^Ew] £7 (1 r 1 £7(2 ) r; - £7 (1 f 1 £7 ( 3 ) £7(if 1 £7(2)7? . (C.18) 

The theory graph contains more lines than constraints, so again, the constraints are compli- 
cated and dependent on the form of the potential. We see as well for the tri-vertex that, due 
to the form of the constraints, the vielbein theory is no longer equivalent to the metric theory 
with the replacement E(ij~ x E{j) —¥ \/ gw -1 gv). It remains to be seen if there exists an equiv- 
alent metric formulation for the vielbein loop-graph theories and theory graphs containing 
higher point vertices, beyond the bi-vertex. 
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